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>-^ . Abstract. The influence of compressibility on the stability of the scaling regimes 

PJ [ of the passive scalar advected by a Gaussian velocity field with finite correlation 

time is investigated by the field theoretic renormalization group within two-loop 
approximation. The influence of compressibility on the scaling regimes is discussed 
as a function of the exponents e and rj, where e characterizes the energy spectrum of 
the velocity field in the inertial range E oc k^^"^^, and rj is related to the correlation 
QQ I time at the wave number k which is scaled as /c"'^+''. The restrictions given by nonzero 

O^l ' compressibility on the regions with stable infrared fixed points which correspond to the 

^— p , stable infrared scaling regimes are discussed in detail. A special attention is paid to the 

f-^ ' case of so-called frozen velocity field, when the velocity correlator is time independent. 

^^3 ' In this case, explicit inequalities which must be fulfilled in the plane e—r] are determined 

^p I within two- loop approximation. The existence of a " critical" value ac of the parameter 

of compressibility a at which one of the two-loop conditions is canceled as a result of the 
competition between compressible and incompressible terms is discussed. Brief general 
analysis of the stability of the scaling regime of the model with finite correlations in 
^ ' time of the velocity field within two-loop approximation is also given. 
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1. Introduction 

One of the main problems in the modern theory of fully developed turbulence is to verify 
the validity of the basic principles of Kolmogorov-Obukhov (KO) phenomenological 
theory and their consequences within the framework of a microscopic model [1, 2]. 
On the other hand, recent experimental, numerical and theoretical studies signify 
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the existence of deviations from the well-known Kolmogorov scaling behavior. The 
scaling behavior of the velocity fluctuations with exponents, which values are different 
from Kolmogorov ones, is known as anomalous and is associated with intermittency 
phenomenon [2]. Even thought the understanding of the intermittency and anomalous 
scaling within the theoretical description of the fluid turbulence on basis of the "flrst 
principles" , i.e., on the stochastic Navier-Stokes equation, still remains an open problem, 
considerable progress has been achieved in the studies of the simplifled model systems 
which share some important properties of the real turbulence. 

The crucial role in these studies is played by models of advected passive scalar fleld 
[3]. Maybe the most known model of this type is a simple model of a passive scalar 
quantity advected by a random Gaussian velocity fleld, white in time and self-similar in 
space, the so-called Kraichnan's rapid-change model [4]. It was shown by both natural 
and numerical experimental investigations that the deviations from the predictions of 
the classical KO phenomenological theory is even more strongly displayed for a passively 
advected scalar fleld than for the velocity fleld itself (see, e.g., [5] and references cited 
therein). At the same time, the problem of passive advection is much more easier to be 
consider from theoretical point of view. There, for the flrst time, the anomalous scaling 
was established on the basis of a microscopic model [6], and corresponding anomalous 
exponents was calculated within controlled approximations (see review [5] and references 
therein) . 

In paper [7] the fleld theoretic renormalization group (RG) and operator-product 
expansion (OPE) were used in the systematic investigation of the rapid-change model. 
It was shown that within the fleld theoretic approach the anomalous scaling is related 
to the very existence of so-called "dangerous" composite operators with negative 
critical dimensions in OPE (see, e.g., [8, 9] for details). Important advantages of 
the RG approach are its universality and calculational efficiency: a regular systematic 
perturbation expansion for the anomalous exponents was constructed, similar to the 
well-known e-expansion in the theory of phase transitions. 

Afterwards, various generalized descendants of the Kraichnan model, namely, 
models with inclusion of large and small scale anisotropy [10], compressibility [11] and 
finite correlation time of the velocity field [12, 13] were studied by the field theoretic 
approach. General conclusion is: the anomalous scaling, which is the most important 
feature of the Kraichnan rapid change model, remains valid for all generalized models. 

In paper [12] the problem of a passive scalar advected by the Gaussian self-similar 
velocity field with finite correlation time [14] was studied by the field theoretic RG 
method. There, the systematic study of the possible scaling regimes and anomalous 
behavior was present at one-loop level. The two-loop corrections to the anomalous 
exponents were obtained in [15]. In paper [13] the infiuence of compressibility on the 
problem studied in [12] was analyzed. In what follows, we shall continue with the 
investigation of this model from the point of view of the infiuence of compressibility 
on the stability of the scaling regimes within two-loop approximation. It can lead to 
sufficient restrictions of the parameter space where the stable fixed points can exist. 
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This, as we shall see rather complicated task, is the first nontrivial step on the way to 
understand the influence of the compressibility of the system on the two-loop corrections 
to anomalous dimensions of the measurable quantities [15]. 

2. Description of the model 

We consider the advection of a passive scalar fleld 6 = 6{x) = 9{t, x) which is described 
by the stochastic equation 

dtO + vAO = uo^e + f\ (1) 

where dt = d/dt, di = d/dxi, z/q is the coefficient of molecular diffusivity (hereafter all 
parameters with a subscript denote bare parameters of unrenormalized theory; see 
below), A = (9^ is the Laplace operator, and /^ = f^{x) is a Gaussian random noise 
with zero mean and correlation function 

{f{x)f{x')) = 5it - t')C(r/L), r = X - x', (2) 

where parentheses (...) hereafter denote average over corresponding statistical ensemble. 
The noise maintains the steady-state of the system but the concrete form of the 
correlator is not essential. The only condition which must be fulfilled by the function 
C{r/L) is that it must decrease rapidly for r = |r| ^ L, where L denotes an integral 
scale related to the stirring. The velocity field v(x) obeys a Gaussian distribution with 
zero mean and correlator [13] 

{v,ix)v,ix')) =Dl{x,x') (3) 

= / ^^ {P^A^) + «Q.,(k)) D\uj, k) exp[-^c.(t - t') + ^k(x - x')], 

where k = |k| is the wave number, u is frequency, d is the dimensionality of the x space. 
In what follows, we shall work with compressible velocity field which is demonstrated 
by the form of the tensor structure of the correlator (3), namely, it consists of two 
parts: the standard transverse projector -Pij(k) = 6ij — kikj/k"^, and the longitudinal 
projector QijCk) = k^kj/k"^ which is related to compressibility. The parameter a > is 
a free parameter. The value a = corresponds to the divergence-free (incompressible) 
advecting velocity field. The function D^ is chosen as follows [12, 13] 

„ ,,3l.4—d—2e—r] 

^^("' ^) = r ^ ,?i • ^ I^y (4) 

The correlator (4) is related to the energy spectrum via the frequency integral 

2 

E{k) ~ k''-' f dcoD^cu, k) ~ ^-^k^-^'. (5) 

J Uq 

It means that the coupling constant Qq (more precisely Qq/uq [13]) and the exponent e 
describe the equal-time velocity correlator or, equivalently, energy spectrum. On the 
other hand, the constant Uq and the second exponent r] are related to the frequency 
ui ~ UQVQk"^'^ which characterizes the mode k. Thus, in our notation, the value 
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e = 4/3 corresponds to the well-known Kolmogorov "five-thirds law" for the spatial 
statistics of velocity field, and rj = 4/3 corresponds to the Kolmogorov frequency. Simple 
dimensional analysis shows that the charges go and mq are related to the characteristic 
ultraviolet (UV) momentum scale A (of the order of inverse Kolmogorov length) by 

^^^^2.+^^ Mo^A". (6) 

In the end of this section, let us briefly discuss two important limits of the considered 
model (3), (4) (see also [12, 13]). First of them is so-called rapid-change model limit 
when Mo ^ oo and Qq = Qq/uI = const, 

D''{u,k)^g',uok-'-''+\ (7) 

and the second one is so-called quenched (time-independent or frozen) velocity field 
limit which is defined by mq — > and Qq = Qo/uq = const, 

D^{lo, k) -^ g'^u^k-''+^-^'7i6{Lo). (8) 

Here the velocity correlator is independent of time in the t representation. 

3. Field Theoretic Formulation of the Model 

The stochastic problem (l)-(3) is equivalent to the field theoretic model of the set of 
fields $ = {9,9',\} (see, e.g., [8, 16]) with action functional 

S{^) = --jdti<f^-,dt2d''^2 t;.(ti,xi)[D];.(ti,xi;t2,X2)]-S(^2,X2) 

+ I dtd'^^e'i-dte-Vidie + uoAe], (9) 

where, in what follows, unimportant term related to the noise (2) is omitted, 6' is an 
auxiliary scalar field, and summations are implied over the vector indices. The second 
line in (9) represent the Martin-Siggia-Rose action for the stochastic problem (1) at 
fixed velocity field v, and the first line describes the Gaussian averaging over v defined 
by the correlator D"" in (3) and (4). 

Standardly, the formulation through the action functional (9) replaces the statistical 
averages of random quantities in the stochastic problem (l)-(3) with equivalent 
functional averages with weight expS'(<l>). Generating functionals of total Green 
functions G(A) and connected Green functions W(A) are then defined by the functional 
integral 

G{A) = e^(^) = fv^ e^(*)+^*, (10) 



where A{x) = {A^, A^ , A"^} represents a set of arbitrary sources for the set of fields $, 
V^ = V6V6'V\ denotes the measure of functional integration, and linear form A$ is 
defined as 

A<^= f dx[A%x)e{x) + A^' {x)e'{x) + A'l{x)vi{x)]. (11) 



Action (9) is given in a form convenient for a realization of the field theoretic 
perturbation analysis with the standard Feynman diagrammatic technique. The matrix 
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{66'), = 
{viVj)a = 



Figure 1. The graphical representation of the propagators of the modeL 



Vj = -ikj = 




Figure 2. The interaction vertex of the model 
(wave-number-frequency representation) . 



of bare propagators is obtained from the quadratic part of the action. The wave-number- 
frequency representation of, in what follows, important propagators are: a) the bare 
propagator {96')q defined as 

{ee\ = {e'e)i = . ] ,, , (12) 

and b) the bare propagator for the velocity field {vv)o given directly by (4), namely 

{v,v,)o = (P,,(k) + aQ,,(k)) D\lu, k). (13) 

Their graphical representation is present in figure 1. 

The triple (interaction) vertex —O'vjdjO = O'vjVjO is present in figure 2, where 
momentum k is flowing into the vertex via the scalar field 9. 

4. Renormalization and RG analysis 

The model under consideration is logarithmic at e = ?7 = (the coupling constants 
go, and mq are dimensionless) , therefore the UV divergences in the correlation functions 
have the form of the poles in e, r/, and their linear combinations. 

The crucial role in the renormalization of the model is played by the total 
canonical dimension of an arbitrary one-particle irreducible correlation (Green) function 
r = ($ ■ ■ ■ $)i_j,,. It plays the role of the formal index of the UV divergence and it is 
given as follows [8, 9] 

dY = d^Y + 2d^ = d + 2-N^d^, (14) 

where A^<j> = {Ng^Ngi^N^} are the numbers of corresponding fields entering into the 
function F, d^ and d^ are the canonical momentum dimension and the canonical 
frequency dimension of the function F, respectively, and summation over all types of 
fields is implied. In what follows, we shall use the definitions of the canonical dimensions 
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of the fields $ as they are given in [12, 13]. It is well-known that superficial UV 
divergences, whose removal requires counterterms, can be presented only in those Green 
functions F for which the total canonical index dr is non-negative integer. 

From the dimensional analysis of the model (see, e.g., [12, 13]), we conclude that 
for any d, superficial UV divergences can exist only in the 1-irreducible functions 
{6'6)i-ir and {6'6v)i_ir. To remove them one needs to include into the action functional 
the counterterm of the form 6'A9 and O'vidiO. Their inclusion is manifested by the 
multiplicative renormalization of the bare parameters go,uo, and uq, and the velocity 
field V in the action functional (9): 

uq = vZy, go = gfi^^'^'iZg, uq = ufi'^'Z^, v -^ Z^w. (15) 

Here, the dimensionless parameters (7,M,and v are the renormalized counterparts of the 
corresponding bare ones, /i is the renormalization mass (a scale setting parameter), and 
Zi = Zi{g, u,a),i = u, g, u, v are renormalization constants. 
The renormalized action functional has the following form 

5($) = --|rftirf'^Xirft2C?''x2t;.(tl,Xi)[D];.(ti,Xi;t2,X2)]-S(^2,X2) (16) 

+ I dt rf'^x e' [-dt9 - Z2Vidie + uZiAe] , 

where the correlator D]^j is written in renormalized parameters (in wave-number- 
frequency representation) 

DUuj, k) = (-P^.(k) + ^g.,(k))^z/V-+u^-'^-^--' ^ ^^^^ 

^■^ ' {iuj + uv^'^k?~^^){—iuj + uv^'^k?~'^) 

By comparison of the renormalized action (16) with definitions of the renormalization 
constants Zj, i = g,u,v (15) we are coming to the relations among them: 

Zu = Zi, Zu = Zi , Zg = Z^Z^ , Z^ = Z^- (18) 

The second and the third relations are consequences of the absence of the 
renormalization of the term with D^ in renormalized action (16). 

The issue of interest are especially multiplicatively renormalizable equal-time two- 
point quantities G(r) (see, e.g., [13]). The example of such quantity are the equal-time 
structure functions 

^„(r) = ([^(t,x)-^(t,x')]") (19) 

in the inertial range, specified by the inequalities /~l/A<<r<<L = \jm (/ is an 
internal length). The infrared (IR) scaling behavior of the function Gir) (for r jl ^ 1 
and any fixed r jL) 

G[r) ~ vfl-'^^[rll)-^^R[rlL) (20) 

is related to the existence of IR stable fixed points of the RG equations (see next section). 
In (20) d^ and dc are corresponding canonical dimensions of the function G, R{r /L) is 
so-called scahng function which cannot be determined by RG equation (see, e.g., [8]), 
and Aq is the critical dimension defined as 

AG = 4 + A.rf^ + 7G- (21) 
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Here 7^ is the fixed point value of tlie anomalous dimension 7^ = fid^ In Zq, where Zq is 
renormalization constant of multiplicatively renormalizable quantity G, i.e., G = ZqG^ 
[13], and A^^ = 2 — 7* is the critical dimension of frequency with '-/^ = 71 which is defined 
further in the text (for more details see, e.g., [12, 13]). 

On the other hand, the small r/L behavior of the scaling function R[r/L) can be 
studied using the Wilson OPE [8]. It shows that, in the limit r/L — > 0, the function 
R{r/L) have the following asymptotic form 

R{r/L)=j:CAr/L){r/L)'^-, (22) 

F 

where Cp are coefficients regular in r/L. In general, the summation is implied over 
certain renormalized composite operators F with critical dimensions Ap. 

In present paper we shall study only the first stage of the RG analysis, namely, 
the infiuence of compressibility of the velocity field on the stability of possible scaling 
regimes of the model. The infiuence of compressibility on the anomalous scaling (the 
second stage of the RG analysis) will be studied in the subsequent paper. 

In what follows we shall work with two-loop approximation. But the calculation 
of higher-order corrections is more difficult in the models with turbulent velocity field 
with finite correlation time than in the cases with (5-correlations in time. It is related to 
the fact that the diagrams for the finite correlated case involve two different dispersion 
laws, namely, uj cc k'^ for the scalar field and uj oc A;^"'' for the velocity field which 
complicates situation even in the one-loop approximation [12, 13]. But, as was discussed 
in [12, 13, 15], this difficulty can be avoided by the calculation of all renormalization 
constants in an arbitrary specific choice of the exponents e and rj that guarantees 
UV finiteness of the Feynman diagrams. From the calculational point of view the 
most suitable choice is to put rj = and leave e arbitrary. Thus, the knowledge 
of the renormalization constants for the special choice r^ = is sufficient to obtain 
all important quantities as the 7-functions, /3-functions, coordinates of fixed points, 
and the critical dimensions. But such possibility is not automatic in general. In the 
model under consideration, it is the consequence of an analysis which shows that in the 
minimal subtraction (MS) scheme, which is used in what follows, all needed anomalous 
dimensions are independent of the exponents e and rj in the two-loop approximation. 
But in the three-loop approximation this dependence can simply appear [15]. 

Now let us continue with renormalization of the model. The relation S{6, 6', v, cq) = 
S^{9,9',\-,e, fi), where cq stands for the complete set of bare parameters and e stands 
for renormalized one, leads to the relation W{A,eQ) = W^{A,e, fi) for the generating 
functional of connected Green functions. By application of the operator "D^ = fidf^ 
at fixed Cq on both sides of the latest equation one obtains the basic RG differential 
equation 

VRGW''{A,e,fi) = 0, (23) 

where I^_rg represents operation 2)„ written in the renormalized variables. Its explicit 
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^H'll 




Figure 3. The one and two-loop diagrams which contribute 
to the self-energy operator Eg'e- 



form is 

I^RG = ^M + 1^9(9, u)dg + (3u{9, u)du - -fu{g, u)V^, (24) 

where we standardly denote V^ = xd^ for any variable x, and the RG functions (the (5 
and 7 functions) are given by well-known definitions and, in our case, using the relations 
(18) for renormalization constants, they have the following form 

7i =V^\nZ„ z = l,2 (25) 

Pg = fxd^g = g{-2e - r/ + 871 - 272), (26) 

I3u = ^^^^,u = u{-r| + '-fl). (27) 

The renormalization constants Zi, and Z2 are determined by the requirement 
that one-particle irreducible Green functions {9'9)i_ir and {9'9\)i-ir must be UV finite 
when are written in renormalized variables. In our case, it means that they have no 
singularities in the limit e, 77 — *> 0. 

The one-particle irreducible Green function {9'9)i_ir is related to the self-energy 
operator S^/e by the Dyson equation 

{9'9)i^ir = -iuj + vqp^ - T.0>0{u,p), (28) 

where the self-energy operator Sg/g is represented by the corresponding one-particle 
irreducible diagrams. In the two loop approximation, it is defined by the diagrams 
which are shown in figure 3. 

On the other hand, the renormalized function {9'9\)i-ir is defined as 

{9'9v,)i-„ = Z2V^ + V,, (29) 

where the function Vj is defined by diagrams of figure 4 (in two- loop approximation). 

Thus, Zi, and Z2 are found from the requirement that the UV divergences are 
canceled in (28), and (29) after substitution z/q = vZy = vZ\. This determines Zi, and 
Z2 up to an UV finite contribution, which are fixed by the choice of the renormalization 
scheme. In the MS scheme all renormalization constants have the form: 1 -(- poles in 
£, ?7 and their linear combinations. As was already mentioned, in our calculations we 
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Figure 4. The one and two- loop diagrams which contribute to the function V^ 



can put rj = 0. This possibility essentially simplifies the evaluations of all quantities 
[12, 13, 15]. The analytical expressions for one-loop diagrams in figure 3 and figure 4 
(in the MS scheme) have the following form 



G 



Sd 



gi^p 



ip 



[1 + u){d - 1 + a) - 2a f fiV' I 



{2nY4u{l + uy 



d 



m 



J 



Gil 



Sd 



m 



a / n 



2e 



(30) 
(31) 



{27tY Au{l + uy d \mj e'' 
where Gip is result for the one- loop diagram in figure 3, and Giy is result for the 
one-loop diagram in figure 4. Here, 5*^ = 27i'^^'^/T{d/2) denotes the d-dimensional 
sphere. The two-loop expressions for the diagrams in figure 3 and figure 4 are rather 
huge, therefore we shall not present their explicit form separately but rather we present 
complete expressions for renormalization constants Zi, and Z2 which have the following 
structure 

Zi = |A, + ^Q5, + a), ^ = 1,2. (32) 
Now using the definition of the anomalous dimensions 71,2 in (25) we obtain 

71 = /i^^lnZi = -2{gA, + 2fC^), (33) 

72 = fJ^d^ In Z2 = -2{gA2 + 2fC2), (34) 

where we denote g = gSd/{27iY. The one-loop contributions Ai and A2 in (33) and (34) 
are defined as follows 

A,= - 
A2 



[l+u){d-l + a) -2a 



a 



d 



Adu{l + u) 



(35) 
(36) 
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and the two-loop contributions Ci and C2 have the form 

Ci = ,.,, ,], , ,, (Cio + «Cn + a^C,2) , (37) 



lQd^u^(l+u, 



3 



where 

((i + 2)(l + m)^ 

_rf-l M(rf-l)(2 + M)(2(M-2)M + rf(2 + 3M)) 

1 + -u a(l + -u)'^ 

(rf - l)(4(-2 + m)m + rf2(-2 + 3^2(2 + u)) + 2d{2 + u{5 - 5u + u"^))) ^^ ,^^, 
+ W^uf ^^' ^^°^ 

^12 = ^Y^^ + ^^Y^(2rfVl + 1.)^ - (« - 3)(« - 1M2 + «) 



+ rf(2 + m)(1 + m(-2 + (m - 6)m))) 

Ci2(l+M) 



+ w^n%il(2(^ - 3)(« - 1)« - ^d'u'il + iz) 



- rf2(_i + u{5 + m(2 + m)(3 + (m - 6)m))) 

+ rf(-2 + m(1 + m(22 + {u- 4)m(2 + u))))) , (41) 

C21 = -{d - l)d{l + uf -{d- l)du{3 + 5u + 2u^)Ho 

2(d-l)d(u-l) , , 

+ (rf - 1)(1 + 2m) (2 + rfM(2 + u))H^ + ^ -f-^ ^/f2, (42) 

C22 = -4rf(l+M) 

+ m(-2(2 + m) + rf(5 - 2m(1 + m) + rf(-l + 2m(1 + u))))Hq 

2(2(1 + m) - rfV(l + Uf + rf(-3 + M + 5^2 + M^) + rf2(l _ y + ^3 ^ ^4)) 

^ d{\\u) ' 

(2 + rf)(l+M) '' ^ ^ 

where we have used the following notation 

d 1 



Hi = 2F1 



1,1-1 + 



2' il + u) 



1 = 0,1,2 (44) 



for the corresponding hypergeometric function 2-^1(0, h]C]z\ = 1 + f^z + \!^(c+i)-i^2 '^^ '^ 
. . .. The functions Bi,i = 1,2 which are introduced in (32) are not important in what 
follows, therefore we shall not define them explicitly. 

5. Fixed points and scaling regimes 

Possible scaling regimes of a renormalizable model are directly given by the IR stable 
fixed points of the corresponding system of RG equations [8, 16]. The coordinates of the 
fixed point of the RG equations are defined by /9-functions, namely, by requirement of 
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their vanishing. In our model the coordinates g^,, -u* of the fixed points are found from 
the system of two equations 

Pg{g*,u*) = Pu{g*,u^) = 0. (45) 

The beta functions f3g and /9„ are defined in (26) and (27). The IR asymptotic behavior 
is governed by the IR stable fixed point which is given by the positive eigenvalues of the 
matrix Q of the first derivatives: 

il, = I 'flf^ "f^lf^ \ . (46) 

V ^^ujdg d(3u/du J 

The infiuence of compressibility on the scaling regimes of the present model in 
one-loop approximation was investigated in [13]. We are interested in the answer on 
the following question: how can the two-loop approximation change the picture of the 
scaling regimes discussed in [13]? 

In what follows, we shall try to study possible scaling regimes in detail. First of 
all, we shall investigate the rapid-change limit: u -^ oo. In this regime, it is necessary 
to make transformation to new variables, namely, w = 1/u, and g' = g/u"^, with the 
corresponding changes in the (3 functions: 

(3g,=g'i-2e + r] + ^,-2^2), (47) 

(3^=w{r]--fi). (48) 

It is well-known that in the rapid change model the higher than one-loop corrections to 

the self-energy operator are equal to zero. On the other hand, the renormalization of the 

velocity field is absent at all as a consequence of the fact that Z2 = 1 aX all orders of the 

perturbation theory. It can be also seen directly by the corresponding manipulations 

with our 7-functions (33) and (34). Therefore, we are coming to the one-loop results of 

[13] (in the rapid-change model limit), namely 

,d — 1 + a , , 

11 = 9 ^ , 72 = 0, (49) 

where again g' = g'Sd/i^nY. 

In this regime we have two fixed points denoted as FPI and FPU in [12, 13]. The 
first of them is trivial one 

FPI: w,=g', = Q, (50) 

with 7^ = 0, and diagonal matrix Vt with eigenvalues (diagonal elements) 

VLi = r], n2=ri-2e. (51) 

Thus, this fixed point is IR stable when 77 > 0, and, at the same time, 77 > 2e. The 

second point is defined as 

2d 
FPU : w, = 0, g: = ——^{2e - r^), (52) 

d — 1 + a 

with exact one loop result 7J = 2e — rj. The corresponding Q matrix is triangular with 
diagonal elements (eigenvalues): 

ni = 2{ri-e), n2 = 2€~ri. (53) 
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It means that this kind of the fixed point is IR stable when rj < 2e together with rj > e. 
The second special case of the present model is so-called "frozen regime" with the 
frozen velocity field. It is obtained from our model in the limit u -^ 0. To consider this 
transition, it is again appropriate to change the variable g to the new variable g" = g/u 
[12]. Then the (3 functions are transform to the following ones: 

/5,» = /(-25 + 27i-272), (54) 

(3u =u{-r] + -fi), (55) 

with unchanged f3 function for parameter u. In this notation, the anomalous dimensions 
7i,2 have the form 

71= ~2i-g"A'; + 2g"'C';), (56) 

72= -2ig"A'^ + 2g"'C'^), (57) 

where, as obvious, g" = g"Sd/{2-nY, and the one-loop contributions are now given as 

K--'-^^. (58) 

and the two-loop contributions C'l and Cg are now 

C'l = Y^ {c'l, + aC'l, + a'C'l,) , (60) 

C'i = ^,[c.C'i, + a'C'^), (61) 



^('0 = ^^^^02 = ^-1, (62) 

C'l, = {d-l)(l- ^!^^Hoi) =l-d, (63) 

C'u = -1 + ^^01 = 0, (64) 



C-2(.-l) (^0.-^^.02)=^, (65) 

C^ = ^(l^iM^^o, + ML±i) ^,, . 4., (66) 

where we denote 

d-2 + 2i 



Hoi — 2F1 



d 



i > 1. (67) 



d-4: + 2i' 

The system of P functions (54) and (55) exhibits two fixed points, denoted as FPIII and 
FPIV in [12]. They are related to the corresponding two scaling regimes. One of them 
is trivial, 

FPIII: u, = g': = 0, (68) 
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with 7i = 72 = 0. The eigenvalues of the corresponding matrix Q, which is diagonal in 
this case, are 

fii = -2€, ^2 = -V- (69) 

Thus, this regime is IR stable only if both parameters e, and rj are negative 
simultaneously. The second, non-trivial, point is 

e 



FPIV: M. = 0, g: 






(70) 



2(A'/ - A'2') 2(A'/ - ^2J 
with exact one-loop relation 7* = 72 + £• After substitution of the corresponding 
quantities one obtains the following expression for the coordinates of the fixed point 



-u* 



0, 9^ 



2de 



1 + 



{d 



1-a 1 + 



d-iy^ ' {d-lf 
The eigenvalues of the matrix Vt (taken at the fixed point) are 

C'i - C'i 



d{d-2) 



2a' 



After corresponding substitutions one has 

e 



Qo 



£-V + l2- 



(^1 = 2£ n 



Qo 



e -1] + 



(d-iy 
ae 



(d-l) 1 



a 



-2a' 



-1+e 



d{d - 2) 
2a\d - 2)d -{d- l)(rf2 - 2) - a{2 + {d - 3)d^) 



(71) 



(72) 



(73) 



(74) 



{d{d - iy{d - 2)) 

The conditions g" > 0, fii > 0, and 1^2 > for the IR stable fixed point lead to the 
restrictions on the values of the parameters e and rj. First, suppose that e < 0. Then 
from the conditions g" > 0, and i7i > one has the following restrictions which must 
be fulfilled simultaneously 

1 + eD <0, 1-eD <0, (75) 

but they cannot be fulfilled at the same time. Thus, our first condition is e > 0. In (75) 
D is given as 



D 



1 



{d 



(rf-1) 1-a 1 + 



2a' 



(76) 



d{d-2)J 

To have g'^ > 0, and r2i > together with e > 0, the following inequalities must be held 
-1<£D < 1, (77) 

which restricts the value of e as a function of the parameter a, and the dimension of 
the space d. In the incompressible case {a = 0) the condition (77) is reduced into the 
simple inequality 

e<d-l. (78) 

In the general case, for each value of d, there exists a "critical" value of a in which 
D = 0. We denote it as ac- In this situation e can be arbitrary, i.e., the condition (77) 
is fulfilled automatically. The value of ac is defined as follows 

2-4d + 3d^-d^ + {4- 16d - 4rf2 + 36d^ - 23ci^ + 2d^ + d^'^ 



ar 



Ad{d - 2) 



(79) 
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For example, for rf = 3 its value is ac = q~ — 0.468. Therefore, in the compressible 
model, the situation is a little bit more complicated as a result of a competition between 
incompressible and compressible terms within two-loop approximation which leads to 
the existence of ac- How does it work? The answer is the following. If we continuously 
increase the value of the parameter a, the region of stability of the fixed point defined by 
the inequalities (77) increases too. This restriction vanishes completely when a reaches 
the "critical" value ac- In this rather specific situation the two-loop influence on the 
region of stability of fixed point defined by condition (77) is exactly zero. Then, if the 
value of parameter a increases further, the condition (77) appears again, and restriction 
on e becomes stronger when a tends, in principle, to infinity. In this limit e ^ 0. On 
the other hand, it must be stressed that in our model only relatively small values of 
a are allowed {a <^ 1). It corresponds to small fluctuations of the density p in the 
system which is supposed in our investigation. In other words, it is supposed that the 
stochastic component of the velocity field of the fluid is much smaller than the velocity 
of the sound in the system (the Mach number Ma <^ 1). 

The last condition on the stability of the IR fixed point is found from the 
requirement to have ^22 > 0. It reads 

,2/ J n\^ (A ^\(A'2 n\ ^Yo i (A '3\a2\ 



as 
rj < € + ■ 



2a\d - 2)d -{d- l)(rf2 _ 2) - a(2 + (rf - 3)^^) 
^^ {d{d-iy{d-2)) 



^0) 



-1 
In the incompressible case it is reduced into the simple condition 

V<e, (81) 

which is held at each order of the perturbation theory. 

In the end, let us consider the most interesting scaling regime with finite value of 
the fixed point for variable u. The coordinates of the fixed point is now defined by the 
requirement of vanishing of the (3 functions which are given in (26) and (27). The fixed 
point value for g = gSd/{2nY is given as 

FPV: 9. = - .^. " .^ - o.?"^^. ^^ (82) 



2(^1-^2) 2(Ai-A 



2 



where the functions Ai, A2, Ci, and C2 are given in (35)-(38), and where the parameter 
u is taken at its fixed point value -u* which is given implicitly by the equation 

-r7 + 7*(M,) = 0. (83) 

Using the exact relations 

7i* = ^, 72 = ^ - ^ (84) 

the expression for the fixed point value of g can be rewritten as a series (expansion) 
of the parameter r^ or a linear combination of r} and e. For example, in [13], where 
the problem was analyzed in one-loop approximation, it was expressed as a function of 
2e — rj (in our notation). In the framework of one-loop approximation it allows one to 
have linear dependence of g^, on the fixed point value of the parameter u. Together with 
another choice of the linear combination of rj and e, namely r] — e it leads also to the 
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simple expression for the fixed point value of u. Thus, the coordinates of the fixed point 
in one-loop approximation are [13] 

2d{l + u^) a r]-2e 

a — 1 + a a— l + arj — e 

It allows, together with the requirement of the positive eigenvalues of the corresponding 
matrix of the first derivatives Q, to find simple conditions for the IR stable fixed point. 
They are defined by inequalities e > 0,e > r], and r] > £ '^~j_~" [13]. 

The situation is essentially more complicated when we are working in two-loop 
approximation. It is given by the fact that now we have nonlinear dependence of g on 
the parameters rj and e, and the expression for the fixed point value of u is now given 
only implicitly in rather complicated expression containing hypergeometric functions. 
Another complication, which defends to analyze the problem in general, is related to 
the fact that contrary to the incompressible case when one has additional condition, 
namely r^ = e, no such condition exists in compressible case under consideration. As a 
result, the analysis of the IR stability of the general case of the present model have to 
be done individually for concrete situation. It is rather cumbersome and it will be done 
in the subsequent work. 

In what follows, let us only give the general analysis of the most interesting case 
when one suppose the relation rj = e. In this situation from the definition of the j3 
functions given in (26) and (27) one obtains the condition 

72* = 0. (86) 

Thus, in this case, the coordinates of the fixed points are given as 

A2{u,) + 2g,C2{u,) = 0, (88) 

but even in this situation the fixed point value of u is defined by complicated implicit 
equation (88) and its exhausted analysis must be discussed separately. 

6. Conclusions 

We have studied the influence of compressibility on the possible IR scaling regimes of 
the model of a passive scalar advected by a Gaussian velocity field with finite time 
correlations by means of the field theoretic RG technique. The possible scaling regimes 
are directly connected to the existence of IR stable fixed points of the RG equations. 
The dependence of the fixed points on the parameter of compressibility and their IR 
stability is discussed. The most attention is paid to the frozen limit of the model where 
inequalities which define the stable IR scaling regimes are found analytically. The 
existence of a " critical" value ac of the parameter of compressibility a at which one of 
the two-loop conditions is canceled as a result of the competition between compressible 
and incompressible terms is discussed in detail. The main conclusion is that for the 
small value of parameter a the region of stability is not restricted considerably. It is 
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also shown that the most general case with finite time correlations of the velocity field 
is more complicated within two-loop approximation and have to be consider separately 
once more. 
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